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ABSTRACT 

In Paper I we showed that steady hne-driven disk wind solutions can exist 
by using "simple" models that mimic the disk environment. Here I extend the 
concepts introduced in Paper I and discuss many details of the analysis of the 
steady/unsteady nature of ID line-driven winds. This work confirms the results 
and conclusions of Paper I, and is thus consistent with the steady nature of the 
ID streamline line-driven disk wind models of Murray and collaborators and the 
2.5D line-driven disk wind models of Pereyra and collaborators. When including 
gas pressures effects, as is routinely done in time-dependent numerical models, I 
find that the spatial dependence of the nozzle function continues to play a key 
role in determining the steady/unsteady nature of supersonic line-driven wind 
solutions. I show here that the existence/no nexistence of local wind solutions 
can be proved through the nozzle function without integrating the equation of 
motion. This work sets a detailed framework with which we will analyze, in a 
following paper, more realistic models than the "simple" models of Paper I. 

Subject headings: accretion, accretion disks — hydrodynamics novae, cata- 
clysmic variables — QSOs: absorption lines 



1. Introduction 

As discussed in Paper I (Pereyra et al. 2004a), accretion disks are commonly believed to 
be present in both cataclysmic variables (CVs) and quasi-stellar objects and active galactic 
nuclei (QSOs/AGN). In both types of objects blue-shifted absorption troughs in UV res- 
onance lines are sometimes present, giving direct observational evidence for an outfiowing 
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wind. Another property that CVs and QSOs have in common is the existence of a per- 
sistent velocity structure in their absorption troughs (when present) over significantly long 
time scales (Paper I). In order for a line-driven disk wind to account for the wide/broad 
resonance line absorption structures observed in many CVs and QSOs, it must be able to 
account for the steady velocity structure that is observed. The 2.5D time-dependent line- 
driven disk wind models of Pereyra and collaborators, both for CVs and QSOs, have steady 
disk wind solutions (Pereyra 1997; Pereyra, Kallman, & Blondin 1997, 2000; Hillier et al. 
2002; Pereyra & Kallman 2003). The earlier ID hne-driven disk wind models of Murray and 
collaborators also find steady disk wind solutions (Murray et al. 1995; Murray & Chiang 
1996; Chiang & Murray 1996; Murray & Chiang 1998). 

However, reports that line-driven disk winds are "intrinsically unsteady" (see Paper I) 
persist, based on the argument that the "unsteadiness" is physically reasonable because of 
the increasing gravity along the streamlines at wind base that is characteristic of disk winds. 

Since the steady nature of CV and QSO wind fiows is an observational constraint, whether 
line-driven disk winds arc steady or not is a significant issue. In Paper I we showed that 
an increase in gravity at wind base does not imply an unsteady wind solution. We also 
developed mathematically "simple" models that mimic the disk environment and we showed 
that line- driven disk winds can be steady. 

In the past, evidence in favor of and against steady line-driven disk wind solutions has 
come from numerically intensive 2.5D disk wind models. However, model difi'crcnccs can be 
numerical in nature. Therefore, a detailed analysis that develops well-defined criteria, and 
applies methods independent of the previous numerically intensive 2.5D models, becomes 
important. The goal of this paper is to extend the concepts of Paper I to estabhsh such 
methods and criteria. In a following paper (Pereyra et al. 2004c, hereafter Paper HI), we 
will apply the criteria developed here to the exact flux distribution of a standard (Shakura 
& Sunyaev 1973) accretion disk. 

I have included gas pressures effects throughout this paper, rather than assuming gas 
pressure to be zero as we did in most of Paper I for reasons of simplicity. It has been argued 
many times in the literature that the assumption of neglecting gas pressure effects (i.e., zero 
gas pressure) in a line-driven wind is reasonable since the main driving force is radiation 
pressure. However, there are four reasons why it may be important to include gas pressure 
effects. 

First, as was discussed by Castor, Abbott, Sz Klein (1975, hereafter CAK75) for the 
stellar case, and as we discuss in Paper I and here, an important property of steady supersonic 
line-driven winds, under the Sobolcv approximation for the line radiation force, is that it 
must have a "critical point." In turn, the position of the critical point determines the exact 
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wind mass loss rates and velocity laws within the model. The critical point is not just a 
mathematical construct that facilitates the calculation of basic wind parameters. As was 
shown by Abbott (1980) for the stellar case, the critical point is the point where the flow 
velocity equals the backward velocity propagation of density perturbations, referred to as 
radiative-acoustic waves or Abbott waves (rather than sound waves). Thus, in a line-driven 
wind the critical point plays a role physically equivalent to the sonic point in a temperature- 
driven wind (Parker 1960). 

However, if one assumes that gas pressure is zero, one finds an infinite family of super- 
sonic solutions of which only one actually has a critical point. The question that naturally 
arises from this specific result is: Does a steady line-driven wind necessarily have to present 
a critical point? The answer is: Yes (for a wind that reaches supersonic speeds under the 
Sobolev approximation for the line force). The argument is relatively simple. As was shown 
by CAK75 for stellar winds, no matter how small the gas pressure may be, the wind solution 
must have a critical point. Therefore, since the zero-gas-pressure case is a limiting case of 
the gas-pressure-included case when gas pressure tends to zero, it follows that only solutions 
with a critical point are "real" or "physical" solutions. Therefore, the requirement of a crit- 
ical point is a gas pressure effect. That is, gas pressure is necessary to identify a unique 
physical solution, which further analysis shows it to be a solution that presents a critical 
point. 

To illustrate the significance of this, if one assumes that gas pressure is zero (rather 
than considering the limiting case when gas pressure tends to zero), the infinite family of 
possible solutions results in wind mass loss rates that may vary from any value arbitrarily 
close to zero up to the wind mass loss rate of the solution which contains the critical point. 
Therefore, the fact that the wind mass loss rate in a supersonic line-driven wind arrives at 
the maximum possible value (that of the critical point solution) is a gas pressure effect ^. 

This means that although the exact value of the gas pressure may not significantly 
affect the actual value of the wind mass loss rate of a steady wind, the existence of gas 
pressure causes the wind mass loss rate to have a unique value that corresponds to the 
critical point type wind solution. Thus, without gas pressure effects, the wind mass loss 
rate would not be uniquely determined, possibly having arbitrarily low values, which in turn 
could have important effects on the evolution of astrophysical systems where line-driven 
winds are present. 



"'^In sections 3.1 and 3.2 of Paper I, where we neglect gas pressure, a unique wind mass loss rate is 
determined because of all the possible solutions only the critical solution is considered. In turn, the reason 
the critical solution is the only physically acceptable one is a gas pressure effect, as can be clearly seen in 
the Appendix of Paper I. 
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A second motive for considering gas pressure effects is that for cases where a steady hne- 
driven disk wind solution is not possible (when assuming gas pressure to be zero), inclusion 
of gas pressure effects may allow a uniquely well-determined solution. An example of this is 
the "S model" of Paper I. If one assumes gas pressure to be zero for the S model, then the 
corresponding nozzle function is monotonically decreasing, and thus the critical point is at 
infinity. Since information, in principle, cannot travel an infinite distance in a finite time, 
a steady physical solution is not realizable for the zero gas pressure case. However, when 
gas pressure effects are included, the corresponding small corrections to the nozzle function 
produce a minimum in the nozzle function at a finite distance (Paper I). Thus, once gas 
pressure effects are considered, the critical point is no longer at infinity, but rather at a finite 
well-defined distance, allowing for a unique physical solution to be found. The fact that the 
inclusion of gas pressure effects may lead a system from not presenting the existence of a 
steady solution to presenting the existence of a steady solution, is a an obvious and more 
than sufficient motivation to include gas pressure effects. 

Third, since evidence in favor and against the existence of steady line-driven disk wind 
solutions comes from numerically intensive 2.5D disk wind models which include gas pressure, 
I am unavoidably led to include gas pressure. If I do not include gas pressure, for example, I 
leave open the possibility that some subtle gas pressure effect at the wind base could generate 
strong fiuctuations down stream causing apparent "intrinsic unsteadiness." 

Fourth, and finally, gas pressure is physically present. For example, in their line-driven 
stellar wind paper, CAK75 did not discuss/justify why they included gas pressure in their 
model, they simply included it since pressure must be present in any hydrodynamic system. 

The specific form of plots, that would further illustrate the general criteria developed 
in this paper, will depend on the specific models being analyzed. However, it is probably 
best to consider a specific model after the general framework is discussed. For this reason, I 
show a series of plots corresponding to the CAK75 model, so as to illustrate and apply the 
concepts and results of this work to a well-known well-studied model. In a following paper 
(Paper HI), we shall apply the concepts discussed here to the case of an exact Shakura- 
Sunyaev disk flux distribution (Shakura & Sunyaev 1973) and compare the results with the 
cataclysmic variable disk wind models of Pereyra (1997) and Pereyra, Kallman, & Blondin 
(2000). 

I present a brief discussion of the steady/unsteady vs. the stable/unstable wind char- 
acteristics in §2. In §3 I present the ID hydrodynamic equations which pertain to the 
line-driven winds models of this work. I extend the definition of the nozzle function given 
in Paper I to include temperature gradients in §4 and discuss the critical point conditions 
in §5. In §6 I develop general criteria for the existence of local steady wind solutions in 
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ID models and discuss the existence of global solutions in §7. In §8 I apply the criteria to 
the well-known and studied CAK75 stellar wind. In §9 I briefly discuss ongoing efforts and 
refer to a following paper (Paper III) in which we will apply the criteria developed here to 
more realistic disk wind systems than those previously analyzed in Paper I. A summary and 
conclusions are presented in §10. 

2. Steady Winds vs. Stable Winds 

In general, a dynamical physical system is said to be steady if and only if there exists 
stationary or time-independent solutions to the equations of motion. A system is said to be 
stable if, in addition to it being steady, arbitrarily small perturbations to the steady solution 
will cither be damped or transported as a wave with constant amplitude by the system. 
Thus, systems with steady solutions that amplify arbitrarily small perturbations are not 
stable. 

Since, in general, physical systems arc in practice subjected to small perturbations, the 
issue of stability is relevant. In particular, if a hydrodynamic flow presents a steady solution 
that is however unstable, then it will be virtually impossible to flnd in nature such a system 
in its steady state. 

Abbott (1980) showed that for the case of ID line-driven wind, under the Sobolev 
approximation for the line force, small perturbations will travel as radio-acoustic waves or 
Abbot waves (rather than sound waves) . The velocity of the radio-acoustic wave is subsonic 
in the direction of the flow and supersonic in the backward direction. Thus, perturbations to 
the steady solution of a ID line-driven will not be amplified but rather will be transported 
through the wind, implying in turn that the wind is stable ^. 

Therefore, for the specific case of a ID hne-driven wind under the Sobolev approxima- 
tion, if a steady solutions exists, then it is also stable ^. 



^In Paper I, as an illustrative alternative to the above argument, and as a consistency check, we showed 
that the steady ID line-driven winds analyzed were stable through numerical simulations. Numerical simu- 
lations, due to the finite precision of the calculations, inherently present small numerical perturbations. 

^Several authors have pointed out that more realistic treatments of line-force that go beyond the Sobolev 
approximation may lead to flow instabilities (e.g., Owocki & Puis 1999). However, the Sobolev approxi- 
mation for the line-force is a standard working assumption in line-driven wind models, which is currently 
being implemented by us and has been implemented in the numerical models that have reported "intrinsic 
unsteadiness" . 



-6- 



3. Equation of Motion 

Following the notation of Paper I, the ID hydrodynamic equation of motion for a sta- 
tionary line-driven flow is: 



1 - 



\ ,dW ^ f A dWy ^^dA db^ 

2W J dz \M dz J dz dz ^ ^ 



which is the equation of motion of Paper I corrected for temperature gradient effects; where 
z is the independent spatial coordinate, W = is the kinetic energy per mass, V is the 

velocity, A is the area that depends on 2;, M = pVA is the wind mass loss rate, p is the 
density, B represents the body forces that corresponds to the gravitational plus continuum 
radiation force per mass, b is the isothermal sound speed, and 7 is the "line opacity weighted 
flux" that also depends on z. 

As in Paper I, I scale the physical parameters by defining a value of Tq, Bq, Aq, and 70 as 
the characteristic distance, gravitational acceleration, area, and line opacity weighted flux, 
respectively. The explicit expressions for the normalized parameters are given in Paper I. 

Introducing the scaling equations, the equation of motion becomes 



/ s\ du ^ f a du\'^ da ds ,. 

(l--)a—^-ga + fa[-—) + 2s— - 2a— , 2 
\ 00/ dx \m dx J dx dx 

which is the normalized equation of motion of Paper I corrected for temperature gradient 
effects; where x is the normalized independent spatial coordinate, uj is the kinetic energy per 
mass, a is the area that depends on x, rh is the wind mass loss rate, g is the gravitational 
plus continuum radiation force per mass, s is the sound speed squared, and / is the "line 
opacity weighted flux" that also depends on x. In equation (2) all physical variables are 
normalized. 

Equivalent to the CAK75 independent spatial variable m, I define the variable g. 



X 

J a{x') 



(3) 



Xo 



where xq is an arbitrary position and qo an arbitrary value. 
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The equation of motion, in terms of variables q and ou, becomes 



A s\duj ^ J. / 1 duj\°' ^ 2s da ^ds 

\ ujJ dq \m dq J a dq dq 

Similar to the CAK75 notation, I define the function h{q) as the sum of the independent 
terms in the equation of motion (i.e., the terms that are independent of cu and duj/dq), 



2s da ds 
a dq dq 



The equation of motion now becomes 



da ds 

ga + 2s-; 2a— 

dx dx 



(5) 



s \ duj , , , ^ / 1 dijj\" 

The problem of existence of a steady solution is thus reduced to determining whether a 
value of m and a normalized function uj{q) exists such that it satisfies the boundary conditions 
and equation (6). A steady solution for the hydrodynamic ID model exists if and only if 
equation (6) is integrable while simultaneously satisfying the boundary conditions. Typically 
the boundary condition is the position of the sonic point. That is, the equation of motion 
must not only be integrable, but additionally there must exist a solution to the equation of 
motion that presents a sonic flow speed at a predefined sonic position (boundary condition) . 



4. Nozzle Function and Critical Point 

To continue the discussion 1 recall the definition of the (3 function given in Paper I, 



(5{uj) ^ 1 - - 

LO 



(7) 
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I extend the definition of the nozzle function to include gradient temperature terms, 



n{q)= q;(1-o;)(^-'^)/"- 



for h{q) < 

(/a)V" 



{ga - [2s/a][da/dq] + 2[(is/dg])(i-")/' 
{ga - 2s[da/dx] + 2a[ds/dx]Y^-'^y'^ 



(8) 



The question of existence of a steady solution is reduced to the question of whether or 
not, upon integration, one can always determine 



dcu dcu . . 



(9) 



Viewing dw/dq as a function of variables q and cu which satisfies equation (6), one can 
divide the q-w plane into five regions depending on whether a solution for dw/dq exists. I 
do this below following the notation of CAK75 for early type stars: 



Region I: uj < s and < : one solution; 



Region II: a; > s and h{q) < and P{oj) fn < n{q) : two solutions; 



Region III: u; > s and h{q) > : one solution; 



(10) 



Region IV: cu > s and h{q) < and ^{00) m > n{q) : no solution; 



Region V: a; < s and h{q) > : no solution. 



As in Paper I, I make the following five assumptions with respect to the solution uj{q): 
uj{q) increases monotonically, the wind starts subsonic, the wind ends supersonic, the wind 
extends towards infinity, and duj/dq is continuous. Following the arguments presented in 
Paper I, and with the additional assumption that asymptotically 
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r a{q) x\q) 

-dx' + qo [i.e., x ^ oo] : < s{q) Soo [> 0] , (11) 
xo ( g{q) a{q) const 

a solution to the ID equations must then have the following sequence in the q — cu plane as 
q increases {x increases): 

Region I: subsonic, h{q) < 0; 

Region II: supersonic, h{q) < , lower branch, /3(a;) m < n{q); (12) 
Region II/Region IV boundary: supersonic, h{q) < 0, critical point, /3((^) rh — n{q); 
Region II: supersonic, h{q) < 0, upper branch, /3(a;) rh < n{q); 
Region III: supersonic, h{q) > . 

I extend the definition of a critical point type solution discussed in Paper I to solutions of 
the equation of motion that satisfy the above conditions [eq. (12)]. Thus, as in the isothermal 
wind case (Paper I) , I have shown for the assumptions given here that a steady solution must 
be a critical point type solution of the equation of motion. 



5. Critical Point Conditions 

As discussed in §4, if a steady solution exists for a supersonic ID line-driven wind 
[with the five assumptions on the u!{q) function described in §4], it must present a critical 
point. The critical point is in the boundary between Region II and Region IV, that is, it is 
supersonic 

U!c> S . (13) 

The sum of the independent terms of the equation of motion [see eqs. (5) and (6)] is negative. 



h{qc) < 



(14) 
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and it satisfies the condition 



(3{uJc) m = n{qc) , (15) 



where the subscript c denotes critical point. 
I define 



and following the notation of CAK75, I define 

F{q, u, u') ^(l-^)co'- h{q) - fa {-\ " {Jf ■ (17) 

Then, the following three conditions must hold at the critical point: 

1) the equation of motion 

F = ; (18) 

2) the critical point at the boundary between Regions II and IV (§4) 

dF , , 

ay"" ■' (1^) 

and 3) the continuity of u' (velocity gradients) ["dF/dq — 0" combined with eq. (19)] 

dF dF 

One thus finds three equations [eqs. (18)-(20)] with four unknowns, namely: qc, cOc, i^'c, and 
the value of rh. 

Therefore, the critical point cannot be uniquely determined by a{q), f{q), g{q), and 
s{q). Thus, when gas pressure effects are included, the nozzle function n{q) cannot by itself 
determine the exact position of the critical point, contrary to the case where gas pressure 
effects are neglected (Paper I). 
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The position of the critical point is determined with an additional model constraint 
which is normally the sonic point position. That is, the equation of motion, upon integration 
from the critical point to lower velocities, must obtain the correct sonic point. 

In §A I analyze in detail the critical point conditions and derive explicit expressions for 
iOc-, uj'c and m as functions of the critical point Qc (eqs. [A15-A17]). Additionally, in §A, I 
also find that in order for to be a critical point (i.e., in order for cUc, <^c' ^o be 

determinable), it must hold that 

h{q^ < and ^ + ^ / , ^ ^ > . 21 

\2sc dq J 1 — a Sc n[qc) dq 

For an isothermal wind {ds/dq — 0) these two conditions for the existence of a critical point, 
reduce to 

> . (22) 

9c 

Although the nozzle function n cannot by itself determine a unique value for the critical 
point, it can constrain the location of the critical point, and in some cases it can be shown 
without additional calculations that a steady solution does not exist (e.g., in an isothermal 
wind with a monotonically decreasing nozzle function). 

The existence of a critical point (i.e., the determination of values for u^., uo'^, and m 
[eqs. (A15)-(A17)] such that the critical point conditions hold [eqs. (18)-(20)]) does not 
imply that the equation of motion [eq. (6)] is locally or globally integrable. That is, the 
existence of a point that satisfies the critical point conditions does not ensure that a local 
(in the vicinity of the point) steady solution exists or that a global (throughout the spatial 
range being considered) steady solution exists. 

In the work of CAK75 for hne-driven stellar winds, the existence of a steady solution 
was proved by finding a solution for the particular case of stellar winds. In the following 
subsection I show that general criteria for the existence/nonexistence of local solutions for 
arbitrary geometries can be established through the nozzle function without integrating the 
equation of motion. 



%, , dn 

c) < and — 
dq 
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6. Requirements for the Existence of Local Steady Solutions 

A local solution in the vicinity of the critical point exists if, upon integration, the solution 
maintains itself in Region II [see eqs. (10) and (12)]. This is equivalent to the condition that 
at the critical point 



P{u!c + Ao;) 771 < 7i{qc + Aq) 



(23) 



where Aq is a variation of q in the vicinity of qc, and Acu is the corresponding variation of 
cu determined upon the integration of the equation of motion. 

Therefore, a local steady solution exists if 



— [P{uj)77i-7i{q)] 



< 



(24) 



I define 



dp ., (fp . d7i „ d'^71 , du d'^uj 

P =1~ ' = ' "'=:r ' "'=TT ' ^ = ^ ' = ■ 

duj dw^ dq dq^ dq dq^ 

(25) 

Prom equation (24), if the critical point conditions hold for a given g, then a local solution 
in the vicinity of that point exists providing 



P"{u) 771 {u'f + P'{u) 7hu" - 7i"{q) < . (26) 

For a given g, the variables cj, cj', and m can be determined through equations (A15), 
(A16), and (A17), respectively. Considering the definition of /3(a;) [eq. (7)], P'{uj) and P"{uj) 
are given by the following equations: 



= ^ ; (27) 

P"{u) = -2-^ . (28) 

The parameter 7i"{q) can be calculated from the functions that define the specific model, 
namely a(g), s(g), g{q)-i and /(g), through equations (5) and (8). Additionally, however. 



the evaluation of equation (26) requires the determination of lo" . This can be obtained as 
follows. A solution to the equation of motion must, of course, be such that the equation of 
motion F = holds [eq. (6); eq. (18)]. Therefore, 



That is 



Thus, in general 



dF , , 

17 = ° ■ 



dF ,dF ,,dF ^ 



dF , JdF 



^"^ ''e, . (31) 

doj' 



But, at the critical point, both the denominator and the numerator of equation (31) are 
equal to zero [eqs. (19) and (20), respectively]. Therefore, 



dF _ ^idF 



< = lim o- (32) 



or 
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Thus, at the critical point 



d 
dq 



'dF' 



d 
dq 



jdF' 



d 

dq 



'dF' 
do? 



= 



(34) 



Therefore, 



+ 



+ 



d^F , d^F 
dudq 



duj \dujdq duP' duj'duj 



d^F 
duj'dq 

|2 



(35) 



d'^F 
du'dq 



d'^F 

du'duj 



+ u'- 



which is a second order equation with respect to oo" . The partial derivates in equation (35) 
can be calculated through equation (17). The values of a;c, o;^, and fn can be calculated 
through equations (A15)-(A17). 

Given ID line-driven wind models with arbitrary geometries, gravitational fields, fiux 
distributions and temperature structures, within the approach of this work, the first step 
in analyzing the existence of steady solutions is to determine for which spatial points the 
critical point conditions as discussed in §5 hold. 

If there are no points where the critical point conditions hold, then a steady solution 
does not exist. On the other hand, if these conditions hold for some of the points, then the 
next step is to determine which subset of points allow a local solution [i.e., satisfy eq. (26)]. 
If there are no points where local solutions of the equation of motion are possible, then a 
steady solution does not exist. 

However, the existence of local solutions does not ensure the existence of a global solu- 
tion. For example, if a global minimum of the nozzle function exists farther out than the set 
of points that allow a local solution [and assuming h{q) < from the set of local solution 
points out to the global minimum], then a global solution does not exist; this is because the 
condition (3{u!) m < n{q) for being in Region II of the q — cu plane (rather than the nonsolu- 
tion Region IV) breaks down before the solution is extended to the nozzle global minimum 
point. At the critical point ^{ujc) Th — n{qc); and (3{u!) is a monotonically increasing function 
[eq. (7)]. 



As indicated above, the local solutions in this work refer to local solutions in the vicinity 
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of a critical point. The requirement of a critical point for the existence of a steady solution 
was discussed in §4. 

Thus, the requirements for the existence of local steady solutions are: 

%e)<0 , (36) 



2sc dqj 1 — a Sc Ti{qc) dq ' 
and 



I3"{ujc) m Kf + I3'{ujc) mJi - v!\q,) < . (38) 



7. Requirements for the Existence of a Global Steady Solution 

The existence of a global solution (a solution that spans throughout the spatial range of 
the model) requires the existence of a critical point (§4), which in turn implies that the critical 
point conditions must hold at the given point (§5), and requires, of course, the existence of 
a local solution (§6) which can be extended from the critical point towards infinity in the 
outward direction and towards the sonic point in the inward direction. 

Additionally, it is required that the global solution be such that upon integration in the 
inward direction the wind reaches sound speed at the sonic point. The exact position of the 
sonic point is a boundary condition of the model (e.g., in the CAK75 stellar wind model the 
position of the sonic point is assumed to be approximately equal to the stellar photospheric 
radius) . 

The set of points for which the critical point conditions hold, and that allow a local 
solution, define the range in which the critical point must be if a global solution exists. The 
exact position of the critical point is determined by initially guessing with a value within 
this range, and iteratively adjusting the critical point position until the correct sonic point 
is achieved when integrating inward. This iterative process to determine the exact critical 
point position in a ID line-driven wind was originally applied by CAK75 in the study of 
stellar winds. Although CAK75 did not present an independent proof of the existence of 
local solutions as 1 do here, in §8 1 show that in the isothermal CAK75 model, for the 
stellar parameters used by CAK75, the points that allow local solutions extend from the 
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photospheric radius to beyond 100 times the photospheric radius. In the original CAK75 
model the iterative process was done over a range of a few integers of the photospheric radius. 

Thus, the requirements for the existence of a global steady solution are (in addition to 
the requirements of existence of a local steady solution) that a critical point exists such that 
upon integration of the equation of motion, these following two conditions hold: 

u;{qs) = s , (39) 

where denotes sonic point, and for the points in the supersonic regions where h{q) < 
(other than the critical point Qc) 

I3{u) m < n{q) . (40) 

8. Application to the Isothermal CAK75 Stellar Wind 

It is well known that when a Sobolev treatment is used for the line radiation force, 
a steady wind solution for early- type stars is obtained (e.g., Owocki & Puis 1999). In 
Paper I we illustrated concepts and notations by applying them to the well-known well- 
studied CAK75 stellar wind, neglecting effects of gas pressure. Here I apply the concepts 
and notations that I present in this second paper, but now including gas pressure effects 
under the assumption of an isothermal wind. 

The purpose of this section in not to shed new hght on the CAK75 model, which is 
already a well-studied model, but rather to illustrate the new concepts introduced in this 
paper by applying them to a problem familiar to the ApJ reader. However, a new result 
is presented, that of the proof of existence of a local solution by analytical means rather 
than by numerical means. This is significant because the critical point, by definition, is in 
the boundary between a solution region and a non-solution region, and therefore numerical 
methods which typically average parameter values in the vicinity of a point to extrapolate 
a function will present difficulties around the critical point and will have to be adapted (at 
the critical point) to the specific case of line-driven winds. 

In other words, subtle numerical methods have to be correctly and consistently intro- 
duced geared specifically to the case of line-driven accretion winds, in order to integrate the 
equation of motion from the critical point. Since the purpose of this series of papers is to 
analyze the steady nature of line-driven disk winds in a form independent of previous nu- 
merical efforts, the fact that I have found an analytical method to determine the existence of 
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local solutions, without integrating the equation of motion, becomes relevant. In particular, 
CAK75 in their paper, did not prove the existence of local solutions, but rather reported 
that they had obtained one through numerical integration. Since the difference between our 
models and those claiming "intrinsic unsteadiness" is a numerical one, how far our analytical 
conclusions can be taken becomes important. 

As in Paper I, I introduce the following characteristic scales for the CAK75 model 



ro = i? 



(41) 



and 



C AttR^ yi^eVthJ 

where R is the photospheric radius, G is the gravitational constant, M is the stellar mass, 
Ke is the Thomson cross section per mass, c is the speed of light, C is the stellar luminosity, 
k and a are the CAK75 line force parameters, Vth is the ion thermal velocity, and F is the 
Eddington ratio given by 



AnGMc ^ ' 

For the CAK75 model the independent spatial variable is the distance r to the center of the 
star, and thus x — r/ R. 

Taking Xq — 1 and qQ — —1, the variable q [eq. (3)] becomes 



q^-- . (44) 

X 

Equation (6) becomes the equation of motion, that is 



s 



duj , , , „ f 1 duj\ 



= + . (45) 



where the function h{q) is given by 



, , , As ^ds 

Hq) = -ga 2 — 

q dq 



(46) 
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and in turn 



ga^l and fa^ ^ . (47) 



Since I am assuming an isothermal wind {ds/dq — 0) 



%) = -l-^ . (48) 

For the isothermal CAK75 stellar wind I implement a set of parameters corresponding 
to an 05f star, namely: 



M = 6OM0 , r = 0.4 , R^13.8Rq , (49) 
and I use the line force parameters also used by CAK75, namely: 



A; = 1/30 and a = 0.7 . (50) 

To study the existence of steady solutions, I first consider the h function. As discussed 
in §5, h must be negative at the critical point. From Figure 1, I find that h is negative from 
the photosphcric height to beyond 100 times the photospheric radius. The sonic radius is 
assumed here to be equal to the photospheric radius. The maximum possible value for the 
critical point position, as constrained by h{q), can be determined by equations (44) and (48). 
However, given the additional assumption that the sonic radius is equal to the photospheric 
radius, I do not expect the critical point to be beyond 100 times the photospheric radius. 

Since I am assuming an isothermal wind, it follows that at the critical point the nozzle 
function must be increasing with position [eq. (22)]. I show in Figure 2 the nozzle function, 
and find that the nozzle function is monotonically increasing from the photospheric height 
to beyond 100 times the photospheric height. Thus, the critical point conditions (§5) hold 
for all spatial points between the photospheric radius and beyond 100 times this radius. In 
other words, for all these spatial points there are well defined, well determined values for cUc, 
cu'^, and m. 

The existence of local solutions in the vicinity of each critical point can be determined 
through equation (26). The expressions P' and P" in equation (26) depend on the critical 
point position through the condition u = oOc and m also depends on the critical point 
[eq.(Al7)]. Figure 3 shows that all the spatial points up to 100 times the photospheric radius 
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present local solutions in their respective vicinities. That is, throughout the aforementioned 
spatial range, the left-hand-side of equation (26) remains negative. 

The global solution is found by determining the exact critical point position by the 
constraint of the sonic point position as discussed in §7, and then integrating the equation 
of motion inwards and outwards. Figure 4 presents the velocity vs. position in the CAK75 
stellar wind for the assumed parameters. As a consistency check for the global solution, in 
Figure 5, I plot functions n and fSrh. In the supersonic region while /i < 0, for a steady 
solution to exist, it must hold that j3m<n for points other than the critical point, and 
j3m = n aX. the critical point. 

By comparing the analysis in Paper 1 of a zero-gas-pressure CAK75 wind with this 
analysis of an isothermal CAK75 wind, I find the following two effects of including gas 
pressure (in addition to unavoidably complicating the analysis). First, in the zero-gas- 
pressure case the position of the critical point is undetermined and all points above the 
photospheric radius become critical points because the nozzle function is constant; however, 
when gas pressure effects are included, corrections to the nozzle function arise and a unique 
critical point is determined through the conditions discussed in §7. Second, in the zero-gas- 
pressure CAK75 wind one has to introduce the additional condition of requiring a critical 
point in order to uniquely determine the wind mass loss rate. 

Without gas pressure effects there are infinite possible solutions, each with a different 
wind mass loss rate. They vary between the value corresponding to the solution that contains 
a critical point, down to wind mass loss rates arbitrarily close to zero. When gas pressure 
effects are included, only one solution exists, which is the one that contains the critical point. 
Thus, when gas pressure effects are included, the physical solution for the equation of motion 
is found without additional conditions other than the equation of motion itself. 

9. Future Application to Line Driven Disk Winds 

This paper is part of an ongoing group effort to model and analyze the line driven 
accretion disk wind scenario for QSOs. Our first 2.5D model results [Hillier et al. (2002)] are 
encouraging in that they are roughly consistent with QSO observational constraints. Among 
the model/observational agreements is that we are finding steady/stable winds within our 
models, and therefore steady wind line profiles. 

However, reports that line-driven disk winds are "intrinsically unsteady" (see Paper I) 
persist, and therefore it becomes relevant to present clear straight forward evidence, indepen- 
dent of previous numerically intensive methods, on whether or not steady/stable line-driven 
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disk wind solutions exist. If the line-driven disk winds were "intrinsically unsteady", then 
one would have to discard the line-driven disk wind scenario on observational grounds. 

After considerable efforts and mathematical analysis of the equations, we have concluded 
that, under the Sobolev approximation, a wind line-driven off a standard accretion disk 
(Shakura & Sunyaev 1973), is steady. 

Given the extensive work that has been done in order to arrive at this conclusion, we 
have decided that it would be best to present it in a systematic form placing emphasis on 
the physical ideas on which the results are based. In Paper I, using simplified models, we 
showed that the increase in gravity at wind base along the streamlines, which is characteristic 
of accretion disk winds, does not imply an unsteady wind. Further we showed in Paper I 
that under reasonable (but simplified) conditions line-driven disk winds can be steady. The 
motivation for using the simplified models was to in turn simplify the mathematical analysis 
and thus focus on the physical principles. 

However, Paper I left one important question unanswered: if one established an equiva- 
lent analysis to more realistic disk models, models such as those currently used to numerically 
simulate 2.5D line-driven disk winds, would there continue to exist steady disk wind solu- 
tions? The answer is: Yes. 

But, in order to derive the answer one obviously has to analyze a more realistic and 
mathematically more complex system. Our approach to presenting the results of the more 
realistic models has been to first establish a more detailed theoretical/mathematical frame- 
work that would in turn allow us to analyze the more realistic models. 

The objective of this paper is to present this theoretical framework, and to illustrate 
it by applying it to the well known CAK75 stellar wind. In a forthcoming paper (Pereyra 
et al. 2004c), we will apply the framework developed to the flux distribution corresponding 
to a standard Shakura & Sunyaev (1973) accretion disk. 

10. Summary and Conclusions 

In Paper I we used "simple" models that mimic the disk environment to show that steady 
wind solutions can exist. Paper I emphasized the underlying physics behind the steady nature 
of line-driven disk winds. The goal of this paper has been to extend the concepts introduced 
in Paper I and discuss important aspects of the analysis of steady /unsteady ID hne-driven 
winds that were mentioned in Paper I, but not discussed in detail. 

Specifically, I show that when including gas pressures effects, the spatial dependence of 
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the nozzle function continues to play a key role in determining the steady /unsteady nature of 
supersonic line-driven wind solutions. The existence/nonexistence of local wind solutions can 
be determined through the nozzle function without integrating the equation of motion, as I 
discuss in detail in §6. This provides a useful numerical test for models aimed at simulating 
line-driven disk winds. 

This work sets a detailed framework with which we will analyze more realistic models 
than the "simple" models of Paper 1. In a following paper (Paper 111), we shall apply the 
framework discussed here to the case of an exact Shakura-Sunyaev disk flux distribution and 
show that if the accretion disk is capable of sustaining the corresponding wind mass flow, the 
wind driven off a standard disk is steady. This is important in that, in turn, it shows that 
the likely scenario for the formation of absorption troughs in CVs is a line-driven disk wind. 
It also shows that a line-driven accretion disk wind continues to be a promising scenario to 
explain the broad absorption lines in QSOs. In the following paper we shall also compare 
the results with the cataclysmic variable disk wind models of Pereyra (1997) and Pereyra, 
Kallman, & Blondin (2000). 

I wish to thank David A. Turnshek, Stanley P. Owocki, Kenneth G. Gayley, Norman 
W. Murray, and D. John Hillier for many useful discussions. This work is supported by the 
National Science Foundation under Grant AST-0071193, and by the National Aeronautics 
and Space Administration under Grant ATP03-0104-0144. 



A. Analysis of the Critical Point Conditions 

The exphcit forms of equations (18)- (20) are respectively: 



and 



s\ , , „ / 1 



a;/ \m' 



(l--)-afa(-Y{ur-'-0 , 



\u! J dq \u! J dq dq \m J 
Prom equations (Al) and (A2), respectively, one then obtains 




Therefore, 

(i-i).'-/. = ^(i-i) 

which leads to 

uj' _ a{-h) 1 
ui 1 — a uj — s 

Equation (A7) can be rewritten in the form 

^ f 1 - = 

a V Lo) 1 — a 
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Comparing equations (A5) and (A8), I have 



or 



ml 1 — a 



fnj fa 1 — a 



(A9) 



(AlO) 



Substituting equation (AlO) into equation (A3), one finds 



1 /a;' 



1 ds f oj' 



2 \uj J 2s dq \uj J 1 — a 2s dq 



a 



-h) d 



In 



(/«) 



l/a 



(_/l)(l-a)/o 







(All) 



Substituting the nozzle function n [eq. (8)] into the natural logarithm of the third term 
of equation (All), one has 



1 /a/ 

2 W 



1 ds fJ 
2s dq \uj 



a {-h) d[\n{n)] 



I- a 2s 



dq 



= 



(A12) 



Assuming that the wind is either isothermal or that the temperature decreases with q 
(i.e., decreases with x), and taking the positive root of the quadratic equation (A12) one 
finds 



1 ds 
2s dq 



1 ds 



a 1 {—h) dn 



2s dqj 1 — a s n dq 
Comparing equations (A7) and (A13) I have 



1/2 



(A13) 



a{—h) 1 1 ds 

1 — a LU — s 2s dq 



1 ds 
2s dq 



a 1 {—h) dn 



1 — a s n dq 



1/2 



(A14) 



and thus 



a{—h) \ 1 ds 
1 — a \ 2s dq 



+ 



1 ds 
2s dq 



+ 



a 1 {—h) dn 
1 — a s n dq 



1/2' 



(A15) 



Substituting equation (A15) into equation (A13), one finds 



s< ^ + 
2s dq 



1 a 1 {—h) dn 

2s dq) 1 — a s n dq 
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+ 



a{—h) 



a 



(A16) 



In turn, substituting equation (A16) into equation (AlO) results in 



1 ds\ a i^h) dn 

2 dq) 1 — a n dq 



1/2 



+ ^ > (A17) 



a 



[cf. eqs. (36), (37), and (39) of CAK75]. 

Therefore, in order for to be a critical point (i.e., in order for uUc, ^'ci ^ to be 
determinable), it must hold that 



/if ~) d ( ^ ds\^ a 1 {—h{qc)) dn ^ ^ 

2sc dq ) 1 — a Sc n{qc) dq 



(A18) 



For an isothermal wind {ds/dq — 0) these two conditions for the existence of a critical point, 
reduce to 



h{q^ < and 



dn 
dq 



> 



(A19) 
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Fig. 1. — The h function for the isothermal CAK75 stellar wind. A necessary requirement 
for the critical point is that the h function must be negative at that point. This figure shows 
that all points ranging from the sonic point [log{x/xs) = 0] to beyond 100 times the sonic 
height [\og{x/xs) = 2] fulfill this condition. The physical parameters used are: M = 60 Mq, 
r = 0.4, R = 13.8i?0 (= photospheric radius), k = 1/30, and a = 0.7 . 
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Fig. 2. — The nozzle function n for the isothermal CAK75 stellar wind presented at two 
different scales. A necessary requirement for the critical point in an isothermal line-driven 
wind is that the nozzle function must be locally increasing (i.e., dn/dq > [dn/dx > 0]) at 
that point. This figure shows that all points ranging from the sonic point [log{x/Xs) = 0] 
to beyond 100 times the sonic height [\og{x/xs) = 2] fulfill this condition. The physical 
parameters used are the same as in Figure 1. 
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Fig. 3. — Local solution existence for the isothermal CAK75 stellar wind presented at four 

different scales: a local solution exists in the vicinity of a critical point x if the expression 
l3" m {uj'y + j3' rhu" — n" (solid curve) is negative. This figure shows that all points ranging 
from the sonic point \\og{x/xs) = 0] to beyond 100 times the sonic height [\og{x/xs) = 2] 
fulfill this condition. In this plot the terms j3" rn{uj'Y (dotted curve), j3' muo" (short dash 
curve), and —n" (long dash curve) are also shown. The physical parameters used are the 
same as in Figure 1. 
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Fig. 4. — Velocity vs. position (radius) for the isothermal CAK75 stellar wind. The critical 
point position is determined with the condition that, upon integration of the equation of 
motion, the correct sonic point position is found. "C" indicates the critical point and "S" 
indicates the sonic point. The physical parameters used are the same as in Figure 1. 
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Fig. 5. — Necessary condition for the global solution existence for the CAK75 stellar wind 
for the critical point shown in Figure 4: upon the integration of the equation of motion, it 
must hold that rh < n{q) at points other than the critical point [when h{q) < and the 
wind is supersonic], and (5{uj)m — n{q) at the critical point. Presented here is the nozzle 
function n (solid curve) and the (5 m function (dotted curve) vs. position at two different 
scales. "C" indicates the critical point. The physical parameters used are the same as in 
Figure 1. 



